Measurement-induced Nonlinearity in Linear Optics 
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PACS numbers: 03.67.-a, 42.50.-p, 03.67.Lx, 03.65.Ta 



I. INTRODUCTION 

In recent years we have the seen signs of a new tech- 
nological revolution in information processing, a revolu- 
tion caused by a paradigm shift to information process- 
ing using the laws of quantum physics [D. Since the 
pioneering work of Feynman Deutsch and Shor 
U a significant effort has occurred worldwide to develop 
the tools necessary to realise such a revolution. There 
are many possible routes and architectures 0, 0| avail- 
able to develop these quantum information processing 
devices. It has long been thought that photons would be 
an extremely strong contender for realising some quan- 
tum information processing circuits • Many of the pho- 
ton's properties, for instance easy manipulation, have 
made them ideal for this. However, for scalable quan- 
tum information processing we require photons to inter- 
act with one another. To achieve such interactions it 
was known that massive reversible nonlinearities would 
be required [j| . Materials giving such large nonlinearities 
were thought to be (and are still) well beyond our ability 
to manufacture. Knill, Laflamme, and Milburn (KLM) 
however found a way to create such nonlinearities using 
only linear optical elements, single photon sources and 
detectors @- More precisely they showed how it is pos- 
sible using such elements to perform conditionally the 
nonlinear transformation 



|*) = c |0)+ Cl |l)+c 2 |2) 
-» co|0)+ Cl |l)-c 2 |2) 



I*') 



(1) 



The optical circuit (depicted in Fig. ^) creating this non- 
linear transformation uses ancilla modes, one prepared 
with a single photon present and the other empty. The 
nonlinearity was induced by definite measurements of the 
presence of the single photon and the vacuum state in 
the appropriate ancilla modes. This insight has reopened 
the door to all-optical quantum information processing. 



Other optical schemes have been proposed alon g th e 
KLM line to generate such sign shifts [H [H ITlTl^|. 
These operations are generally conditional in nature. By 
this we mean the transformation only works when the ap- 
propriate measurement results are obtained at the ancilla 
detectors. While this would seem to limit the viability of 
the information processing, it is straightforward however 
by using a teleportation-based protocol to turn such non- 
deterministic operations into deterministic ones 0, Il5j| . 
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FIG. 1: Schematic setup of the KLM circuit for generating 
a nonlinear sign shift using 3 beam-splitters, a single photon 
source and single photon resolving detectors. 



There have been a number of key experiments demon- 
strating elements of linear optical information processing 
[Trj Il7l lis|. These have generally focused on the tech- 
nology necessary to perform single-qubit rotations and 
CNOT gates. Such gates are well known to be sufficient 
to perform universal computation (they are the minimum 
set required). From these primitive elements interest- 
ing devices such as Quantum Repeaters yj| and single- 
photon quantum non-demolition detectors ^| can be 
created. In this paper we wish to shift the focus slightly. 
Instead of using only these primitive gates we will investi- 
gate what operations can be constructed from linear ele- 
ments, single photon sources, and detectors. This shift is 
analogous to the shift in classical computing from a RISC 
(reduced instruction set) architecture to the CISC (com- 
plex instruction set) architecture. The RISC-based archi- 
tecture in quantum computing terms could be thought of 
as a device built only from the minimum set of gates 
while the CISC-based machine would be built from a 
much larger set, a natural set of gates allowed by the 
fundamental resources. 
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Our primary focus in this paper will be on the opera- 
tions that can be constructed from the linear optics set. 
We show how to construct general operators that can be 
applied to the required input states. We further indicate 
what operations are easily constructed and what are po- 
tentially difficult, illustrating our constructive procedure 
with examples from one- mode and two- mode situations. 
Our constructive procedure can easily be applied to mul- 
tiple modes. The inputs to the computational modes do 
not need to be restricted to qubits only: the operations 
can be applied onto qudits and continuous variables just 
as easily. 

This paper is organised as follows. In Sec. [H] we will 
derive some general expressions necessary for the con- 
struction of useful nonlinear operations. In Sec. II I II we 
will be concerned with single-mode operations, followed 
by two-mode operations in Sec. IIVI Until then, we as- 
sume perfect beam splitters and detections which is an 
oversimplification, indeed. We will therefore focus on the 
effects of absorption and non-unit detection efficiencies 
in Sec. I VII before drawing some conclusions in Sec. IVllI 
Some useful formulae regarding permanents of unitary 
matrices can be found in the Appendix. 



II. GENERAL BEAM-SPLITTER 
TRANSFORMATION 

In order to introduce the notation we will be using 
throughout the paper we will briefly review the most ba- 
sic features of quantum-state transformation by a lossless 
beam splitter. We refer the reader to the extensive litera- 
ture for details [2(j ■ Every (lossless) beam splitter can be 
thought of as a unitary operator on the level of photonic 
creation and annihilation operators of the incoming and 
outgoing fields, i.e. 



Wall = Ad . d 



b 



(2) 



where U is a unitary operator and A the associated uni- 
tary matrix [A € SU(2)]. The transformation matrix A 
consists of the transmission and reflection coefficients T 
and R and can be given in the form 



A = 



T R 

-R* T* 



(3) 



Unitarity of A requires |T| 2 + |i?| 2 = 1 which leads to the 
usual definition of the beam splitter 'angle' if by writing 
\T\ = cos</2, \R\ = sin ip. The unitary operator U can be 
given in several equivalent forms, two of which are the 
following: 



U 
U 



id ' #d 



e 

jmi e - 



A = e' 19 



(4) 
(5) 



The effect of the beam splitter cannot only be described 
by transforming the photonic operators, but equivalently 



by transforming the quantum state g as 

gout = U Qi a U^ . 



(6) 



Noting that the input density operator g; n can be writ- 
ten as a functional of photonic creation and annihilation 
operators, Q[ n — Qm[a, df], the quantum-state transfor- 
mation can be represented as 



Pout — Qh 



g- m [A+d,A T dt] , (7) 



that is, the state transforms with the inverse operator 
pH I22] ] . On the level of quantum states we thus have to 
perform the replacements 



a 

At 



A+d, 
A T dt. 



(8) 
(9) 



We will use Eq. @ extensively throughout the paper. 

Suppose we were given an input state with N modes 
with the associated creation and annihilation operators 
labelled by ai , i — I...N. Additionally, we have a 



N 



supply of M auxiliary modes labelled by j 
1, . . . , N+M. Then, a general unitary transformation on 
all the modes maps a* h-> A T d t , A G SU (N + M). What 
we mean precisely by SU (N+M) is a unitary operator on 
the level of photonic creation and annihilation operators 
in N + M dimensions. In what follows, we will only make 
use of the unitarity of the corresponding matrices and 
will not further elaborate on the actual underlying group 
structure. In order to construct our quantum operations 
we will use the decomposition of an arbitrary element of 
the group SU(N) into at most N(N - l)/2 U(2) group 
elements, i.e. beam splitters \2?i . 

First, let us define our notation. By 10)®^ we mean 
the tensor product state |0)i |0) 2 ■ • ■ |0)jv- Let the input 
state now be given in a functional form as 



liM = /(a},...,a] v )|o>® JV 

and the auxiliary state in product form as 



N+M 



IV'aux) 



(4) 

j=N+l V m 3- 



(10) 



(11) 



Here nij is a non-negative integer that represents the 
number of photons initially in the mode j. Finally, the 
state we project on shall be denoted by 



proj/ 




(12) 



where rij represents the number of photons in the pro- 
jected mode j. The output state after mixing at the beam 
splitter network and projecting onto IV'proj) looks then as 
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;u ix / 
JV+M 



M<g>, 



n +l v^i^i 



n< /N+M \ ™ l J /N+M N+M \ 

( E A«alJ / ( E AafiJ, E AwaJ ) |0>«"+* 



(13) 



\ fc=i 



/=i 



What we see here is that the effect of the beam splitter network is to generate the desired mixing of the photonic 
creation operators of signal and auxiliary modes. Now we make use of the ordering formula well-known from bosonic 
operator algebras (see, e.g., 0,113) 



d 



[a,F(a,a^)] = —F(a,a^) 



to rewrite the output state as 

N+M 



|Vw) oc 



(oi n 




N+M \ 

, E AlN& l 



\oy 



iN+M 



(14) 



(15) 



Furthermore, we expand the function f(a\, . 
Taylor series as 



N 



f{a\, ■ ■ ■ , ojv) = E 



(a!)'" 



Pl,...,PN 



= 1 



. , a' N ) in a 
(16) 



■"-pi,...,p N 

Cm ', 



is constrained in such a way that 
v | 2 = 1. In that way we obtain the 



where c 

N 

E 

pi,...,PN = l 

action of a SU(7V + M)-network in a quite general way. 
In general, this can be a laborious task. In order to see 
the structure behind it, let us focus first onto single- mode 
signal states. That is, the input state will be 



= M)|o) 



(17) 
(18) 



and the network will represent an element of the group 
SU(iV + l). 

In what follows we will restrict ourselves to the im- 
portant special case when our resources consist of single 
photons and single-photon detectors. In this case, we can 
derive a number of interesting results. Let us first start 
with a very simple (and in fact well-known) example, a 
single beam splitter. Feeding a single photon in one input 
arm of the beam splitter and measuring a single photon 
leaving one output port of the beam splitter, we have in 
fact created the conditional non-unitary operator [using 
Eq. ©] [U 



Y 



(h\U\l 2 ) =T^- 1 [|T| 2 - n x \R\ 2 ] (19) 

acting on some signal state |^i n ) (see Fig. [2}. This is a 
very special result and probably the simplest non-unitary 
operator one can actually generate. This conditional op- 
erator has already been realised in an experiment |27| 
where it is called 'quantum-optical catalysis'. 



IVin> 
ll> 



SU(2) 

V J 



l ( Kout> 

ll> 



FIG. 2: Schematic setup for generating the simplest non- 
unitary conditional operator with a single-photon input and 
a single-photon detection. 



In the following we will present some results on the 
general structure of conditional non- unitary operators. 

Proposition 1: Let us suppose all N auxiliary modes 
are prepared in single-photon states, and all TV detectors 
measure vacuum. This is equivalent to acting with an op- 

rator ~ I a[ j on the signal state (left figure in Fig. [3J| . 



Proof: The auxiliary and detected states are 



N+l 



0I\ 



l^dct) - \0)' 



(20) 



i=2 



The conditional (un-normalised) output state is therefore 

ivw) oc e ■ c " 




/N+l 

E A «4 

\fc=i 



5 (4) 



A 




(21) 
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Apart from normalisation (or success probability) , which 
depends on the chosen input state, the output state is 
proportional to the TV-fold application of the creation 
operator. rj 

In complete analogy, we can prove the following: 

Proposition 2: Let us suppose all N auxiliary modes 
are prepared in the vacuum state and each of the N detec- 
tors measures a single photon. Then, this is equivalent to 
acting with on the input state (right figure in Fig. [2J • 



l¥in> 

n> 



n> 



SU(N+1) 



IV„ut> 

I0> 



I0> 



I0> 



I0> 



IVout> 

SU(N+1)|_^[V ll> 



ll> 



FIG. 3: Adding (subtracting) photons to (from) the signal 
mode by subtracting (adding) the corresponding number of 
photons from (to) the auxiliary modes. 



Proof : Again, let us first write down the auxiliary and 
the detected state: 



N+l 



KW> = |0) 8JV , iVdet)=n«ii 



Acting on the input state gives 



(22) 



|VWt) oc ^2 



N(g>/ 



(N+l \ /N+l \ "* 

/N+l „ \ /N+l \ m 

Ei^(n I A..)Arr»(4)"-*',o, 

(23) 



'N+l \ 

[] A 4l A* 1 of |Vi„> , 

V i=1 J 



where in the last line we have repeatedly made use of the 
formula 



(at) p |0) = ^-a(at)^ |0) 



1 



which immediately follows from the commutation rela- 
tions of the photonic operators. This proves that, in- 
deed, measuring N photons from an TV-mode auxiliary 
vacuum input is equivalent to acting N times with the 
annihilation operator on the signal state. rj 

Propositions 1 and 2 show how to generate arbitrary 
powers of creation and annihilation operators. In fact, 
one could have already guessed the general form of these 
operators by recalling that the network is represented 
by an element of the compact group SU(A + 1). Com- 
pactness of the group translates into photon-number 
conservation which is why adding (subtracting) N pho- 
tons from the auxiliary modes must end up as subtract- 
ing (adding) photons from (to) the signal mode. Note 
that in both cases only the matrix elements An or Ah 
(i — 2, . . . , N+ 1), respectively, appear. This means that 
the network decouples into a sequence of N disconnected 
beam splitters. That is already the minimal number of 
beam splitters necessary for the generation of the wanted 
operators. 

The next step consists of showing how powers of the 
number operator can be realised. In fact, an obvious 
way would be to combine the results from Propositions 1 
and 2 and to construct an alternating network producing 
sufficient numbers of creation and annihilation operators. 
This might not be the most sensible way to do. In fact, as 
we will see later, the following result has much stronger 
implications for the construction of interesting quantum 
operations. 

Proposition 3: Measuring single photons in all N de- 
tectors from a supply of N single-photon auxiliary state 
amounts to multiplying the input state with a polynomial 
of iVth degree in the number operator, P/v(^i) (Fig. EJ. 
Proof: We will only sketch this proof and calculate the 
highest power of h\ and leave the remaining terms for an 
interested reader to calculate. Given that we choose the 
auxiliary and detected states of the form 



JV+l 



l^aux) = n«ll°>^' 
i=1 
N+l 

l^dct) = l[a]\0)® N , 



(25) 



k=2 



(24) the output state can be written in the following way: 



oc 



E 



<o| 



/N+l q 



N+l /N+l 
j=2 \i=l 



'N+l \ 

E A ™ i6 « ) 

vn=l / 
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+1 \ / N+1 

n a « n a - 

j=2 I \n=2 




hi—N 



11* 



'N+1 

3=1 ier 



(26) 



In the first term the factorial hi ! / (ni — N) ! is a polyno- 
mial of order N in hi and thus the desired result. All 
other terms (not written except for the last, in lowest 
order in hi) contain lower-degree polynomials psj . This 
proves the assertion. rj 



IVin> 

n> 



n> 



SU(N+1) 



l¥out> 
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FIG. 4: Generating polynomials of photon-number operators 
by single-photon inputs and detections. 

The simplest example of this proposition is a single 
beam splitter the result of which we have already seen in 
Eq. J5J). However, with the above propositions, we can 
immediately generalise our considerations to obtain the 
following results: 

1. Given that the following for ancilla and detected 
modes 



|V'a 



1 



, ®N-\-Af 



IV'det) - |1) 8JV ®|0)® M ' 

the output state will be 

IVw) « (a\) M P N (hi)\^ in ). 

We immediately see that this procedure has allowed 
us to act on the input state with the creation op- 
erator (a\) M . 

2. Analogously, with 

IVw) = \1)® N ®\0)® M 

l^dct) = \1)® N + M 

the output state will be 

IVw) cx P N {hi){ai) M \1>in) ■ 

In both situations we have, with the aid of linear optics, 
single photon sources and detectors, been able to operate 
on the input state \tp[ n ) with both af 1 and (a\) M . Let 
us now turn our attention to single-mode operations that 
are of interest in connection with quantum information 
processing. 



III. SINGLE-MODE OPERATIONS 

From now on we will focus onto the generation of uni- 
tary operators which are of utmost importance for most 
quantum information processing tasks. For all unitary 
operators it is easy to define the success probability, since 
unitary operators leave the norm of a quantum state un- 
changed. Since these operators Y are prepared condi- 
tionally, the success probability is just 



P: 



success 



(27) 



for any (normalised) state vector \ip). 

We can derive some interesting results about these uni- 
tary operators. For example, let us suppose our input 
state is a single-mode state consisting only of elements 
in the zeroth and first Fock layer. It is clear that all 
operations on |V>in) of the type 

|Vin) =co|0)+ci|l) ^co|0) + e^ci|l) (28) 

can be realised with a probability of p = 1, since uni- 
tary operations simply consist of phase shifts of the |1) 
state. A special example with ip = tt is the Pauli-<7 Z . 
Going one step further we may ask what the conditions 
are for generation of unitary operations on single-mode 
states with up to two photons. It is reasonable to as- 
sume that we would need at least an SU(3)-network, that 
is, two auxiliary modes. In fact, we find that every uni- 
tary single-mode operator acting on states with up to two 
photons, separately in each Fock layer, can be generated 
by an SU(3)-network with two single-photon inputs and 
two single-photon detections. In order to show that, let 
us first calculate the conditional operator for the Sub- 
network with l^aux) = iV'dct) = 1 1 1) - We get 

FlVin) = per A(l|l)|0)+per A|l) 
+ (2A llP er A - A^per A(l|l) + 2A 12 A 21 A 13 A 31 ) |2) . 

(29) 

It is known that the range of per A (as a function of all 
its relevant parameter) is the unit disk in the complex 
plane [2t| (see Appendix). In fact, so is the range of any 
principal sub-permanent per A(i\i). This can be seen 
from the decomposition of an SU(3)-matrix in terms of 
a product of three SU(2)-matrices [23| which themselves 
have a range spanning the unit disk. Therefore, it is 
immediately clear that we can again generate any phase 
e zipi between the states |0) and |1). As for the two-photon 
Fock layer, we can rewrite the coefficient in Eq. I|29|l to 
obtain a condition on the matrix A as 

per A(l|l) [e^ + A 2 U - 2A 11 e^ 1 ] = 2A 12 A 21 A 13 A 3 i , 

(30) 
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where e* V2 is the phase shift between |0) and |2). The 
modulus of the rhs of Eq. (|3U[) can be shown to be 
bounded from above by 8/(27|An| 2 ) by noting that 
Y[i Aii is the product of the elements of a unit vector. 
Noting also that the principal sub-permanent per A( 1 1 1) 
can take any value across the unit disk we can conclude 
that Eq. (|30|1 has always a solution. This in turn means 
that every unitary single-mode operator acting within 
Fock layers on states with up to two photons can be 
generated by an SU(3)-network with two single-photon 
inputs and two single-photon detections which was to be 
proven. The probability of success is |per A(l|l)| 2 . It is 
also possible however to create certain phase shifts with 
the necessity for two ancilla photons. For instance, in Q 
it was shown that a sign shift on the |2) Fock state only 
is possible with the ancilla state 1 10) . 



IV. TWO-MODE OPERATIONS 

In order to do something useful in terms of quantum 
information processing, we have to operate on two modes 
simultaneously. This can be done in more than one way. 
For example, one can simply generalise the theory pre- 
sented above for a single signal mode to more than one 
signal mode. It turns out that this is not a very trans- 
parent way. We will follow another route instead and de- 
compose the two-mode operation into three subsequent 
steps: 

1. combine the two modes at a beam splitter, 

2. act on both modes separately, 

3. and recombine the modes at another beam splitter. 

The effect of the beam splitters is to mix the modes and 
to make them accessible for a single-mode operation in 
such a way that we can apply the result in Sec. IIIII 

A. The 'controlled-phase' gate 

We will illustrate this statement with an example. 
Consider the two-mode operator C v acting on qubits. 
Its truth table is 



|00) -> |00), 

|01) - |01), 

|io> - |io), 

111) e iV |ll) 



(31) 



In terms of photon creation and annihilation operators 
the operator C v can be represented as 



C v = 1 - (1 - e iV )nin 2 . 



(32) 



Now let us assume that we mix the signal modes at a 
symmetric beam splitter. The operator C v acts only in 



the two-photon Fock layer. Then it is very easy to see 
that with (nonlinear) single-mode operators N% = 1 — 
|(1 — e l¥, )n«(nj — 1), i = 1, 2, we achieve a transformation 
of an input state 

\M = c 00 |00) + c O i|01) + cio|10) + c n |ll) (33) 

into 

C^in) = coolOO) + coi|01) + cio|10) + c n e iv \ll) . (34) 
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FIG. 5: Controlled-phase gate with single-photon detectors 
only. 

The nonlinear operator needed on both modes are 
polynomials of second degree in the number operators hi 
and can thus be prepared conditionally with two auxil- 
iary modes prepared in single-photon Fock states on each 
side followed by double single-photon detection. Hence, 
the overall requirements are four single-photon sources, 
eight beam splitters, and four single-photon detectors. 
The generic network is shown in Fig. [3] The detectors 
all measure single photons. We can write down the con- 
ditional operator as 

Y\i> in ) = per A(l|l)c |0) + per A Cl |l) 

+ [2A 12 A 21 A 1 3A 31 + 2per A - A 2 lP er A(l|l)] c 2 |2) . 

(35) 

The success probability is |per A(l|l)| 2 . Numerically, we 
find values up to p S ucccss ~ 0.24 in each interferometer 
arm. 

However, it turns out that there is an even simpler net- 
work with only six beam splitters and two single-photon 
sources |l2|. It has the disadvantage, though, that one 
needs two vacuum detectors which are hard to make (and 
which are pretty inefficient). The corresponding network 
is shown in Fig. EI The set of beam splitters fed with vac- 
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FIG. 6: Controlled-iTz gate with single-photon and vacuum 
detectors. 
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uum states act as conditional phase shifts. In summary, 
we find that the beam splitters must satisfy 



arg^ji) = -argTj 0> , (36) 
\T\i)\ = 0.476, (37) 
|T| 0> | = 0.87, (38) 

which gives a success probability of p succoss ~ 0.23 in each 



arm, hence a total success probability of w 0.05. 

Let us remark that the controlled-(7 2 investigated by 
Ralph et al. 01 f ans the same category as that de- 
scribed in our Fig. [5] The difference is that one of the 
single photons in each arm of the interferometer is re- 
placed by the vacuum state and the single-photon detec- 
tor by a vacuum detector j^3| , respectively. This network 
corresponds to the following conditional operator: 

Y\& m ) = A 22 c |0)+per A(3|3) Cl |l) 

+ (2A 12 A 21 A 11 +A 22 A 2 1 )c 2 |2). (39) 

The probability of success is |A 22 | 2 . One needs to satisfy 
the set of conditions 

P erA(3|3) = A 22 , (40) 
2A12A21AU + A 22 Af l = -A 22 , (41) 

from which it immediately follows that An = 1 — y2. The 
maximal value |A 22 | 2 can take under the constraints 140|l 
is then indeed 0.25 which is why the gate in Ref. [IJ is 
indeed optimal. 

B. The 'swap' gate 

A somewhat more interesting operator is the swap op- 
erator S in the sense that here we encounter the first 
example of an operator that needs fewer resources than 
one would expect when considering CNOT and single- 
qubit rotations as building blocks for quantum circuits. 
It is known that it can be made from three CNOT oper- 
ators (equivalent to controllcd-<7 z gates with attached 
Hadamard gates). Acting on qubits, one can write the 
photonic-operator version of it as 

S = fiin 2 + (ni — l)(n 2 — 1) — a|a 2 (ni — 1) — a 2 ai(n 2 -l) . 

(42) 

Let us see how the single-mode version of S can be de- 
rived. It is immediately clear that we have to act on 
the single-photon Fock layer only. It turns out that the 
nonlinear single-mode operators are 

Ni = l + 2ni(ni-2), (43) 
A> 2 = 1, (44) 

which means that we do nothing on mode 2, and we 
act with a polynomial of second degree in hi on mode 1. 
Therefore, we would need only two single-photon sources, 
four beam splitters, and two single-photon detectors. 



However, the operator N%, when acting on Fock states 
\n), is nothing but a single-mode phase shift (— l)™ 1 . 
That is, the whole network collapses into a single tt- 
phase plate in one arm of the Mach-Zchndcr interfer- 
ometer leaving us with just two beam splitters and one 
phase plate. This gate is remarkable in the sense that it 
is also unconditional, that is, it works deterministically 
with unit probability which makes it rather special. 

These two simple examples show a general principle 
of constructing these networks. Both operators have in 
common that they act only within a specific Fock layer 
(S: one photon; C v : two photons). One then projects 
out all those Fock layers which are not affected by the 
operator. This leads to the polynomials in the number 
operators. The design of the polynomial coefficients in 
each case depend on the specific operation one wants to 
achieve. 



C. General considerations 

A general conclusion can already be drawn from the 
results on one- and two-qubit operators: It is highly de- 
sirable to rewrite the quantum information network in 
such a way that the actual computation can be made as 
long as possible in the same Fock layers. Every crossing 
to another layer (cf. the Pauli operators a x and a y ) re- 
quires additional resources which might not be necessary. 
This leads us to state our main result of this paper: 

Theorem: The generic operations that can be done eas- 
ily and effectively with linear optics are operations within 
the same Fock layers. Let M be the number of signal 
modes we want to operate on. Any M-qubit gate act- 
ing within Fock layers can be constructed with the help 
of generalised Mach-Zehnder interferometers with M in- 
put and output ports (2M-ports for short) and at most 
M conditional operators generating polynomials in the 
number operator of at most Mth order (equivalent to 
SU(Af + l)-networks). 

Proof: The proof of this assertion is now straightforward. 
Any operator acting within Fock layers can be written as 
a polynomial of at most Mth order in all photon number 
operators. The 2M-port mixes all the M input modes in 
such a way that we are left with a tensor product of M 
operators inbetween the 2M-ports conditionally generat- 
ing polynomials of at most Mth order in the individual 
photon number operators. rj 

This result shows how to construct these operations 
in an algorithmic fashion. That is what we mean with 
'easy'. Since there is no inherent exponential scaling of 
the success probability with respect to the number of 
modes (qubits) we act on, there is a good reason to call 
them also 'effective'. 

Unfortunately, not all two-qubit gates can be written 
in terms of a Mach-Zehnder interferometer and appro- 
priate single-mode operations. Perhaps the most notori- 
ous example is the CNOT gate. Although similar to the 
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controlled-<7 z , there is no way to find an interferometric 
setup that 'disentangles' the two modes in such a way 
that there existed single-mode operators that performed 
the sought task. The proof of this statement goes along 
the following lines: Let us call U(<p) the beam splitter 
operator that rotates the qubit axes by an angle ip [see 
Eq. (J5J; a Mach-Zehnder interferometer would consist of 
a succession of two of these operators with opposite an- 
gles]. Here, we seek a transformation of the following 
type: 

IVw) = E%)(iVi ® N 2 )U(cp'M n ) := (45) 

with the two (conditional) nonlinear operators N\ and 
N2. A lengthy but straightforward calculation shows that 
the operator sandwiched between the beam splitters does 
not have tensor-product structure and thus cannot be 
regarded as single-mode operators. In order to show that, 
we use a matrix technique. Let us define a basis vector 
|e) as 



(|00),|10),|01),|11),|20),|02)). 



(46) 



Then, the input state (V'in) can be written as 
iV'in) = c hJ e )- I n this basis the vector = 
(cod) cio, coi, en, 0, 0) transforms as 



cout -U(^)(N 1 ^N 2 )UMc i! 



(47) 



where the matrices U(y) etc. are the matrices corre- 
sponding to the operators U(ip) etc. in the basis |e) 
(these are not to be confused with the beam splitter or 
transformation matrices used earlier on). For example, a 
beam splitter is represented in this basis by the matrix 



Ufa) 



(I 









T R 

-R* T* 












-V2R*T V2RT 
T R 2 

(48) 

with \T\ = cos ip and \R\ — s'm(p. The tensor product of 
the two single-mode operators looks in this basis like 



\T\ 2 -\R\ 2 

V2RT 
-V2R*T* 



J 



Ni<g)N 2 = 



/ (#1)00(^2)00 (#1)01(^2)00 
(iVi)io(JV 2 )oo (^1)11(^2)00 
(#1)00(^2)10 (^1)01(^2)10 



V 



(49) 



/ 



It is then relatively straightforward to show that there 
exists no solution to Eq. (|47|l with a matrix of the 
form (|49|l that produces an output vector c^ ut = 
(coo,cio,cii,c O i,0,0). 

Therefore, in order to build a CNOT gate, we would 
have either to refine our approach to include more gen- 
eral interferometric setups (for which the original Knill- 
Laflamme-Milburn proposal is an example) or sandwich 
a controlled-^ gate between two Hadamard gates which 
we will show in the next section to be rather expensive. 



V. CROSSING FOCK LAYERS 

Equipped with the knowledge about generating anni- 
hilation and creation operators, we can start working on 
realisations of other operations that are harder to do but 
nevertheless needed to construct general quantum net- 
works. By our Theorem, the 'easy' operations are those 
that act within the same Fock layers. It is much harder 
to find suitable networks for operators that enable us to 
cross Fock layers 0]. The obvious choice consists of 
looking at single-qubit rotations first, i.e. the represen- 
tations of the Pauli operators in the Fock basis, 



|0)(1| + |1)(0|, 

-(|o)(i|-|i)(o|). 



(50) 
(51) 



The construction of the corresponding photonic opera- 
tors is almost obvious, once one takes care of the fact 
that one must not leave the Hilbert space of the qubits. 
Then it is clear that we have to choose 



a 
1, 



a^n-l), 
+ a t (n - 1)] 



(52) 
(53) 



In order to proceed further, we need a well-known result 
from quantum-state engineering. 

Proposition 4: Suppose one wants to generate the 
quantum state 



fc=0 



(54) 



Then one needs n single-photon sources, at most n 
coherent-state sources, and at most In beam splitters 
and detectors. 

Proof: The proof of this proposition follows closely the 
result in (sj where it has been shown that the state \ip n ) 
can be generated by successive single-photon additions 
and coherent shifts. The trick is to rewrite the state as 

n 

hu=n( at - a fc)i°> ( 55 ) 

k=l 

which is nothing but a decomposition of the polynomial 
in into its root factors, where the a£ are the roots of 
the polynomial. ■— ] 

Having generated the state \ip n ), one can go ahead and 
imprint it onto another state by mixing at a beam split- 
ter. That leads neatly to 

Proposition 4a: The polynomial 



k=0 



(56) 



can be made to act upon a signal state by mixing the 
state V n \0) and the signal state at a single beam splitter. 



9 



Proof: Let us assume that the signal state is again of the 
form 



l<U=£4=(alrio) 



(57) 



Mixing jV'in) and "P«|0) at a beam splitter, conditional on 
the second output being found in the vacuum state, we 
obtain after a short calculation 



\1>o 



(58) 



fc=0 



from which we see that the coefficients have to be suffi- 
ciently rescaled to achieve the desired goal. rj 

In the same manner one can generate polynomials of 
annihilation operators by projecting onto an engineered 
state. Combining both processes opens up the opportu- 
nity to generate arbitrary polynomials of creation and an- 
nihilation operators. However, this might not be the best 
choice since doing quantum-state engineering of higher- 
order polynomials is, as we have seen, an expensive task. 
Therefore, it might be advantageous to circumvent the 
problem of leaving the Fock layers of zero and one photon 
by projecting back onto this subspace after performing a 
simplified version of the desired quantum operation. For 
this, we introduce the KILL operator K as 



K = 1 



1 



h(n — 1) 



(59) 



which, being a second-order polynomial in the number 
operator, requires two single-photon sources, two beam 
splitters and two detectors. The Pauli operators can then 
be written as 



<t v = K—(a — a^) . 

i 



(60) 
(61) 



With the theory presented above, we could go ahead 
and generate superposition states |0) + |1) with the help 
of Proposition 4a, superpose them onto the signal mode 
and perform a projection measurement onto a similar 
state. However, we will present a slightly different and 
more elegant method of achieving this purpose. Instead 
of preparing two copies of the superposition of vacuum 
and a single photon, we could prepare a Bell-type state 
~ |0,0) + A|l, 1) by the following method. Let us take a 
two-mode squeezed vacuum state of the form 



|TMSV) = yjl- q 2 ^q n \n,n) 



(62) 



ra=0 



and perform a Procrustean [32, |33| entanglement con- 
centration by acting on one mode of it with a first-order 
polynomial of the number operator as explained in the 
example JjJ. For appropriately chosen transmission co- 
efficient T of the beam splitter, we can generate in the 
limit q — > the state 



|$(A)> 



|0,0> + A|1,1>] 



(63) 



to arbitrary accuracy in the trace-norm and for arbitrar- 
ily chosen A (details of this procedure can be found in 
|34|)- Using this state as the auxiliary-state source in 
an SU(3)-network that projects onto |1,0), we derive the 
following operation after applying the KILL operator: 



c |0) + ci|l) -> A 2 ici|0) + A per A(3|l)c |l) 



(64) 



Choosing | A21 1 = |A per A(3|l)| with an appropriate 
phase relation immediately leads to the desired Pauli op- 
erators. 

A remark concerning the usage of continuous-variable 
states as resource is of order here. In the described ver- 
sion of the Pauli operators we inject a two- mode squeezed 
vacuum state into our network. This seems a simple and 
elegant method for getting the desired result. In fact, we 
cannot see a way around the usage of continuous- variable 
states at all, since even for the creation of the superpo- 
sition |0) + |1), by Proposition 4, a coherent-state source 
is needed to displace the photon creation operator a) . A 
similar conclusion was reached by Lund and Ralph 

Another very important single-qubit operation is the 
Hadamard gate, defined by 

10) - -±=(|0} + |1)), (65) 

|1> - ^(|Q>-|1». (66) 

This can also be written in operator form as 

F=^(|0) + (-lf|l>) , (67) 

where the number operator is the one from the signal 
state! That is, we swap signal and auxiliary states in 
the sense that we first produce a superposition of |0) 
and |1) and act conditionally on it with the signal state. 
Effectively, the Hadamard gate becomes a (controlled) 
a z operation on the (auxiliary) superposition state (|0) + 
|1))/V2. In fact, one can rewrite the operator H as 



H 



1 



v/2 



(|0> + (l-2nin 2 )|l)) 



(68) 



which is effectively a two-mode operator. This is pre- 
cisely the controlled-<7 z where we the second output is 
left unmeasured (sometimes called the DUMP 'gate'). 
However, leaving something unmeasured usually means 
to trace over the possible outcomes which will destroy 
the purity and coherence of our desired operation. The 
way around this problem is to act on the resulting signal- 
mode output with an operator 1 + (which can be pre- 
pared according to Proposition 4a) and then to project 
onto the single-photon Fock state. 

Form this rather complicated construction we observe 
that the Hadamard gate, and consequently also its multi- 
mode extension, the quantum Fourier transform, are the 
hardest of all gates under investigation so far. This result 
impacts the generation of gates that actually make use of 
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similar layer-crossings as the CNOT gate. For these type 
of operations it seems that the constructive algorithm we 
have presented in this article is not immediately applica- 
ble and this problem requires further investigation. 



VI. LOSSY BEAM SPLITTERS AND 
NON-PERFECT DETECTORS 

So far, we have restricted ourselves to perfect linear 
optics, i.e. non-absorbing beam splitters and detectors 
with unit efficiency. In practise, to achieve this situation 
is a hopeless task. Instead, we have to make do with 
absorbing linear optical elements and non-perfect detec- 
tors. What this amounts to in terms of constructing our 
gates will be described in this section. 



A. Kraus decomposition 

We derive the Kraus decomposition of a lossy beam 
splitter. It is known that an absorbing beam splitter 
represents a unitary evolution in the extended Hilbert 
space of field and device modes. The unitary operator 
can be written as 133 



where we have used the definitions 



U = exp 
where we use the notation 



a = 



a 



(69) 



(70) 



Assume now the device to be initially in its vacuum state 
|03,04). Then we can write the density operator of the 
output field as 

£ u \ = Tr(°) [0 (^ ) |03,0 4 )(03,0 4 |) L>t] (71) 

and evaluate the trace in the coherent-state basis as 



(F) 

Pout = £2 / dfasdfaiE^^Qll'E 1 



1 



(F) A t 



(72) 



where we have defined the Kraus operators E a3j0l4 as 

E a3 , ai = (a 3 , a 4 |£7|0 3 , 4 ) . (73) 
They can be further simplified by using the relation (3|J 

in . 

(74) 



71=0 



n 



by writing 

(a3,a4|^|03,0 4 ) = (a 3 ,a4|exp 



-i [a 



a 



|0 3 ,0 4 ) 



: [a* (A- l)al" :, 
03, a 4 V -L-^ , ' J 3 ,04 

n=0 



Tl=0 



[at(T-l)a-a+SC" 1 Ta] r 



-e 



e -«6 t * T o e -a+SC- 1 Ta e -|a+( 



(75) 





( T 


A = 


\ — SC 


C = 


Vtt+, 


S = 


Vaa + , 


T = 


e -7* T ; 


g|a 3 ,a 4 ) = 


a a 3 , a 4 ) 



A 



e~**, (76) 



(77) 

(78) 
(79) 
(80) 



Therefore, we obtain the result that the Kraus operators 
for the absorbing beam splitter are 



E n 



g-io 1 * r o e -a + SC~ 1 T& - 



11) 



We can easily check that these operators become unitary 
when absorption can be disregarded as T becomes uni- 
tary (and therefore 3?t hermitian), and S vanishes. The 
integration over (03, a 4 ) can then be performed and gives 
unity. What we also see is that these Kraus operators in- 
deed correspond to an absorption process for which the 
factor exp[— a+SC _1 Td] is responsible. 



B. Non-perfect detectors 

Second, we model a non-unit detector efficiency r\ by 
replacing the projector |n)(n| by an appropriate POVM 

|n)(n| - ft(n) = £ ( fc V (1 - . (82) 



This method does not take care of possible dark counts 
but reflects the fact that direct photon counting may give 
values for the photon number n that actually came from 
higher Fock states \k), k > n. This POVM is sometimes 
modelled by a perfect detector preceded by a beam split- 
ter with appropriately chosen transmissivity |T| 2 = 77. 



C. Example: a single beam splitter 

Let us consider a somewhat artificial example which 
nevertheless shows what happens when absorption 
and/or non-perfect detectors are present. Suppose we 
were to implement the Pauli-<7 Z gate with a single beam 
splitter, a single-photon source and a single-photon de- 
tector (note that this could have been done deterministi- 
cally with a phase plate) . We start off with a signal mode 
in a state co|0) + ci 1 1) and mix it with a single photon. 
The effect of the absorbing beam splitter is to produce a 
mixed state that can be written in the form 



-(F) 
fi'out 



|^ in (T))(^„(T)| + |0(A))(0(A)| (83) 



where \ipi n (T)) is the state transformed with the (non- 
unitary) transmission matrix T and \cf>(A)) is a contribu- 
tion that solely comes from the absorption matrix A. We 
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do not give the rather lengthy expression here. Instead, 
we immediately give the result for the non-normalised 
density matrix after applying the POVM (|8*2")l as 



f?out,l 



??IVw)(Vw| 

+4r ? (l- ?7 )| Cl | 2 |T 12 | 2 |T 22 | 2 |0)<0| 



I 2 | 

-f7|ci| 2 (|T 22 Mn 



T12M21 



|T 22 M 12 +T 12 M 22 | 2 )|0)(0| 



(84) 



with the wanted output state 



|Vout) = c T 22 |0) + ci(T u T 22 + T 12 T 21 )\1) (85) 

and the matrix M = SC X T. Eq. JSlJl has three parts: 
The first line is the wanted outcome in which the trans- 
mission matrix can be chosen to give the desired answer. 
The second line comes from the inefficient detector, hence 
the POVM introduced in Eq. JS^Jl, whereas the last two 
lines are the contributions due to the lossy beam split- 
ter, reflected in the appearance of the matrix M that 
contains the absorption matrix. The last expression can 
be simplified using the fact that MM+ = 1 - TT + to 
obtain 



£out,i = f7|Vw)(Vw| 

+4r 1 (l- V )\c 1 \ 2 \T 1 
+v\ci\ 2 [\T 22 \ 2 + \T 12 \ 2 -4\T 12 \ 2 \T X 



21 ir 22 |lo)(o| 



-|T n T 22 + T 12 T 21 1 2 ] |0)(0|. (86) 

This expression shows that it is only necessary to know 
the experimentally accessible transmission and reflection 
coefficients of the beam splitter that make up the ma- 
trix T. Now we make use of the fact that we actually 
wanted to generate a Pauli-cr z gate, meaning that we set 
in Eq. JESJ T X1 T 22 + T 12 T 21 = -T 22 . With that we finally 
obtain for the (still unnormalised) output density matrix 

&mt,l = ?7|T 22 | 2 <J z |?/>i n )(V'in|<5' z 

+4 ?7 (l~r,)| Cl | 2 |r 12 | 2 |T 22 | 2 |0)(0| 
+V\ci\ 2 [|T 12 | 2 -4|T 12 | 2 |T 22 | 2 -3|T 22 | 

The success probability for perfect operation is p, 
|T 22 | 2 . A note of caution is appropriate here. Since we 
have fixed T 22 already, by reciprocity we have also fixed 
Tn = T 22 = T. For single-slab beam splitters that fixes 
T\ 2 = T 2 i = R, too, so that we are left with essentially a 
single number determining the fidelity of our desired gate 
operation. To be more precise, note that |T| 2 + \R\ 2 + 
\A\ 2 = 1 (setting \A\ 2 = |A n | 2 + |A 12 | 2 = |A 21 | 2 + |A 22 | 2 ), 
and suppose that T G R. Then we immediately have that 
R 2 € R, and choosing argi? = tt/2 we arrive at 



(87) 



success 



T 



(88) 



With this choice for T 22 = T we finally get 
&ut,i = r?(2 - \A\ 2 - y / 3^2\Aj)a z \^ in )(^ in \a z 



+ri(l V )\ci\ 2 (\A\ 4 - 3 + 2V3-2L4P) |0> <0| 

+?7 | C1 | 2 |A| 2 (1-L4| 2 )|0)(0| (89) 

which now only depends on two parameters, the absorp- 
tion coefficient \A\ of the beam splitter and the detector 
efficiency 77. Again, the first line is the desired result, the 
second is due to the non-perfect detector, and the last 
line is the contribution of the absorption. Two special 
cases are notable here: 

1. Without absorption (\A\ = 0), the third line in 
Eq. (|89|l vanishes and the numerical coefficient in 
the second line takes the value of 2y/3 — 3 « 0.464. 

2. With perfect detectors (rj = 1), the second line van- 
ishes and we are left with a contribution |A| 2 (1 — 
\A\ 2 ) to the vacuum from the last line. 

In principle, one could define a (state-dependent) gate 
fidelity or use some more elaborate definition such as an 
average fidelity integrated over all possible input states 
(with respect to some Haar measure) but this is beyond 
the scope of this article. 



VII. CONCLUSIONS 

In this paper we have shown a constructive mechanism 
for generating arbitrary operators using only linear op- 
tics, single-photon sources, and single-photon detectors. 
We have focused our attention primarily on one-mode 
and two-mode situations, though the approach is easily 
extended to multimode situations. We have shown what 
operations are easy and what are potentially difficult. 
Operations that cause a change in the Fock layers (for 
instance the Hadamard operator) are generally difficult 
but not impossible. While the generation of the operators 
is generally conditional on certain measurement results 
in the ancilla modes, the operators can be made deter- 
ministic using various teleportation protocols. Finally 
we hope this paper shows the power in building the re- 
quired operations from the fundamental resources rather 
than fundamental gates. The SWAP operation illustrates 
this point extremely well. From fundamental gates, three 
CNOTs are required to build such an operation, however 
from fundamental resources, only two beam splitters and 
a phase shifter are necessary. This approach open a new 
way to think about operation generation. 
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APPENDIX A: PERMANENTS OF UNITARY 
MATRICES 

Here we recall some elementary properties of perma- 
nents, mainly taken from the only available monograph 
on this subject |2^. The permanent of an (nxn)-matrix 
A is a generalised matrix function, defined as 

n 

perA= ]T Y[A^ (Al) 

WeS™ 1=1 

where S n is the symmetric group of cyclic permutations. 
Note that the determinant of a matrix is similarly de- 
fined with the only difference of a factor of (—1) appear- 
ing in all terms depending on the character (even or odd) 
of the permutation. The permanent of a matrix gener- 
ically appears in counting problems, i.e. combinatorics 
and graph theory. In our case it is the probability am- 
plitude of detecting the state after an input state 
of the exactly the same form has been transformed by 
an SU(7V)-network. In that sense, it naturally appears 
here as well since the combinatorial problem is here to 
(re-)distribute iV single photons among N single-photon 
detectors. 

The Marcus-Newman theorem states that the follow- 
ing inequality holds for all (mxn)-matrices A and (nxm)- 
matrices B: 

|per AB| 2 < per AA*per BB* . (A2) 



An immediate consequence is that (setting B = l), if U 
is unitary, then 

|per U| < 1 . (A3) 

Note that this condition also follows immediately from 
the probabilistic interpretation given above. Eq. (|A3|) 
tells us that the range of the permanent of a unitary 
matrix lies in the unit disk in the complex plane. In fact, 
the same conclusion can be drawn for the permanents 
of principal submatrices of unitary matrices by recalling 
that a unitary matrix consists of rows (or columns) of 
orthogonal unit vectors. For example, let us consider 
per A(l|l) of Ae SU(3). We have 

|per A(l|l)| = |A 22 A 3 3 + A 23 A 32 | . (A4) 

Since |A 23 | < y/1 = |A 22 | 2 and |A 32 | < y/l - |A 33 | 2 , we 
know that 



|perA(l|l)| < |A 22 A 33 | + |v/(l-|A 22 |2)(l-|A 33 |2)| 
= | cos 9? cos 8| + | sin<^sin0| 
= |cos(<p±0)| < 1. (A5) 

Similar relations hold for per A(2|2) and per A(3|3) and 
indeed for all permanents of submatrices of unitary ma- 
trices. 
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